^-dependence of QCD at Finite Isospin Density 

Max A. Metlitski^B and Ariel R. Zhitnitsky 1 '! 

Department of Physics and Astronomy, 
University of British Columbia, Vancouver, BC, Canada, V6T 1Z1 
(Dated: February 2, 2008) 

We probe the ^-dependence of QCD at finite isospin chemical potential [ii using 
the effective chiral Lagrangian approach. The phase diagram in the 9, \xi plane is 
constructed and described in detail in terms of chiral and pion condensates. The 
physics at 6 ~ it is investigated in both the normal and superfluid phase. Finally, 
the behaviour of the gluon condensate at finite [ij is computed. 

PACS numbers: 



I. INTRODUCTION 

The ^-parameter of gauge theories has long attracted attention as it is a probe of the 
topological properties of the theory. In almost every context, from pure Yang- Mills theories 
to QCD the 9 dependence of the theory is highly non-trivial and, frequently, non-analytic. 
In particular, in QCD with two flavours and equal non-zero quark masses it is believed 
that the so called Dashen's phenomenon [J, 0, E3, 0, HI El - a first order phase transition, 
characterized by spontaneous breaking of CP, occurs at 9 = ir. 

In this paper, we investigate the influence of finite isospin chemical potential fii on the 9 
dependence of two flavor QCD. Besides pure academic interest, the main physical motiva- 
tion for such a study is the attempt to understand the cosmological phase transition when 6, 
being non-zero and large at the very beginning of the phase transition, slowly relaxes to zero, 
as the axion resolution of the strong CP problem suggests. Of course, in real world, we are 
mostly interested in the effects of 9 on matter at finite baryon, rather than isospin, density. 
Indeed, if isospin asymmetric matter presently exists in nature (say in neutron stars), it is 
accompanied by a large baryon density. However, analytical control over QCD is absent at 
moderate baryon density and appears only at asymptotically large baryon chemical poten- 
tial, where one expects the color-superconducting state to be realized [7f. Nevertheless, one 
may resort to QCD-like theories, such as N r = 2 QCD at finite baryon densitvP.l9l.llo| and 
N c — 3 QCD at finite isospin density [ill. E3 |. where analytical control is present, to gain 
some insight into real dense QCD. 

Due to the axial anomaly, the 9 parameter of QCD is intimately tied to the quark 
mass matrix and may be incorporated into the effective chiral Lagrangian^, HI- There 
also exists a well-known procedure for including the effects of finite /i/ into the QCD chiral 
Lagrangian [III Il2^ . We, thus, expect that we may adequately describe QCD at finite isospin 
density and 9 ^ in the effective Lagrangian approach, as long as Hi is much smaller than 
the mass of the lightest non-Goldstone boson (in QCD, the mass of the p meson, m p ). 
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Using the above approach, we obtain a wide range of information about the phase diagram 
of two flavor QCD in the ///, 9 plane. We show that the transition to the superfluid, isospin 
breaking, phase occurs at /zj equal to the 9 dependent pion mass, m n (9). This implies that 
for fixed /zj of order of the pion mass, the 9 dependence of the theory becomes non-analytic. 
Two second order phase transitions, accompanied by a jump in the topological susceptibility, 
occur as 9 relaxes from 2n to 0. 

We compute the 9 dependence of chiral and pion condensates, as well as (iGG) and the 
topological susceptibility, in normal and superfluid phases. We find that the 9 dependence in 
the superfluid phase near 9 = it is much smoother than in the normal phase. In particular, 
for m u = ma, we show that the first order phase transition across 9 = ir present in the 
normal phase, disappears in the superfluid phase. 

Finally, we discuss a few 9 unrelated issues. Most importantly, we compute the de- 
pendence of the gluon condensate {^^G^G^ 11 ) on the isospin chemical potential in the 
superfluid phase. The gluon condensate decreases with density near the normal to superfluid 
phase transition, but, counter-intuitively, increases for m n <C /!/ <C m p . We also evaluate 
a novel vacuum expectation value, which appears in the superfluid phase: (iu^olad). This 
density, being nonzero even at 9 = 0, nonetheless has never been discussed in the literature 
previously. This density, itself, breaks the isospin symmetry, and so may be considered as 
an additional order parameter. 

We note that the above agenda has also recently been implemented to study the properties 
of N c = Nf = 2 QCD in the presence of non-zero 9 at finite baryon and isospin density. 
Most of the results of the present study are in direct correspondence with the work 13[ • This 
is a consequence of the fact that the chiral Lagrangians describing N c = 3, Nf = 2 QCD and 
the pion sector of the N c = 2, Nf = 2 QCD are identical. Besides adapting the work|l3[ to 
the N c = 3 context, we presently discuss in some detail the theoretically interesting case of 
exactly degenerate quark masses, which was not analyzed in [13j. 

We hope that the results of this study would be of interest for lattice simulations. Indeed, 
the determinant of the Dirac operator is real and positive in QCD at non-zero isospin 
chemical potential and 9 = 0. The determinant remains real at /i/ ^ 0, 9 = n. Thus, we 
hope that the /i/ dependence of the gluon condensate and the topological susceptibility at 
9 = 0, can be explicitly checked on the lattice. This is a unique chance to study the gluon 
degrees of freedom and their dependence on light quark masses. The corresponding study 
might be important for the extrapolation procedure which has to be used in order to achieve 
the chiral limit. Moreover, we hope that the disappearance of Dashen's phenomenon at 
9 = Ti in the superfluid phase can also be confirmed by lattice simulations. 

II. THE CHIRAL LAGRANGIAN 

The low energy dynamics of Nf = 2 QCD are governed by the chiral Lagrangian for the 
pion field U G SU(2). A well known procedure exists to incorporate into this Lagrangian 
the effects of a finite 9 parameter [2jj3j. A method for introducing a finite isospin chemical 
potential \ij is also well-developed 12j . To lowest order in quark mass and derivatives, 
the chiral Lagrangian reads, 



C = ^Tr^f/V^) - ZReTr(MU) (1) 
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where the flavor covariant derivatives are defined as, 

V U = d U - ^ 7 [r 3 , U], ViU = diU (2) 
V f/t = doU+ l -^[U\r\ = dirf (3) 

We work in Euclidean space. Here the 9 parameter of QCD has been incorporated directly 
into the quark mass matrix, 

e/N f ( m u 



M = e~ l °l N f '7 " (4) 

We keep m u ^ m d on purpose: as is known m u = m d is a very singular limit when one 
discusses 9 dependence, see below. The coefficient S is determined by the chiral condensate 
in the limit m — > + , 9 = 0, fii = 0, 

- <^>o (5) 



2N f 

as will be confirmed below. In our notations the chiral condensate includes the sum over all 
flavors, (tjji/j) = J2 f 

Due to pseudo-reality of SU (Nf = 2), one may, to this order in chiral perturbation theory, 
incorporate all effects of 9, m u , m d into a common real quark mass via a redefinition, 

U = LUR\ L = Ftf = e iar3/2 (6) 
(m u + m d ) cos(6 l /2) 



cos a 



v / (m M + m d ) 2 cos 2 (#/2) + (m u - m d ) 2 sm {6/2) 



(m u -m d )sm(9/2) 
sin a = (7) 

y/(m u + m d y cos 2 (fl/2) + (m u - m d ) 2 sin 2 (#/2) 

Our parameter a is related to the commonly used Witten's variables <p u , cj) d [2], via, 

<P U = 9/2 -a, 4> d = 9/2 + a (8) 
(fiu + 4>d = 9, m u sin <p u = m d sin 4> d (9) 

After such a transformation, the Lagrangian takes the form, 

£ = IfiTriVjrVjP) - m(9)E ReTr(U) (10) 

with, 

m{9) = i ({m u + m d ) 2 cos 2 (#/2) + (m u - m d ) 2 sin 2 (#/2)) ~ 2 (11) 



III. PHASE DIAGRAM 



Our next step is to find the classical minimum of the effective Lagrangian (fTU|) to deter- 
mine the phase diagram. First let's study the theory at zero chemical potential and fixed 9. 
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The classical minimum, is then given by, U = 1, and the lowest lying excitations correspond 
to a triplet of pions, with 9 dependent mass, 

ml(e) = Mffmm (12) 

J IT 

The pion mass acquires a dependence on 9 through the effective quark mass parameter 
m{9) ([IT]). As we shall see, the whole phase diagram turns out to be determined by the 
parameter m v {9). We note that m n (9) reaches its maximum at 9 = and minimum at 
9 = 7i. Moreover, for m u = m^, 9 = tt, vanishes to first order in m q . 

Now let's turn on finite /!/. At fixed 9, the phase diagram contains two phases: normal 
and superfluid. The transition from the normal phase to the superfluid phase occurs at the 
critical chemical potential [ij = m n (9). In the normal phase, (U) = 1. In the superfluid 
phase, the U(l)i symmetry is spontaneously broken and, 



(U) = A(0) + iy/l - A(^) 2 (r 1 cos + r 2 sin 0) (13) 

where the variable labels the U(l)i degeneracy of the vacuum, and we have introduced 
the parameter A to describe both the normal and superfluid phase, 

{1 normal phase 

=5J& superfluid phase. < 14 » 
Pi 

As expected, at 9 = we reproduce the known results 11, At 9 ^ the phase diagram 



looks the same as at 9 = 0, with the important replacement, m 2 — > rn 2 K {9). This is a very 
natural conclusion. Indeed, at 9 ^ pions still carry isospin number. Hence, their energy 
is lowered at finite isospin chemical potential. As soon as \xi reaches the vacuum pion mass 
m n (9), Bose-condensation occurs leading to spontaneous breaking of U(l)i symmetry. 

Quantitatively, the 9 dependence of the Goldstone mass m- K {9) implies that the transition 
to superfluid phase is shifted to a smaller chemical potential compared to 9 = 0. In the 
limiting case, when m u = and 9 = tt, the transition occurs in the vicinity of ji = (see 
Section IV for a more precise discussion). For physical values, = 7MeV, m u = 4MeV, 

the transition at 9 = tt occurs at /i = ( j^Zl ) * ~ 70 MeV. 

We now wish to describe the phase diagram in terms of different condensates and densi- 
ties. This can be achieved by the standard procedure of introducing sources into the chiral 
Lagrangian. We find that chiral condensates depend on fij, 9 in the following way, 

1 9 1 9 

(uu) = -{#) A(fi) cos(- - a), (dd) = -(^} X(9) cos(- + a) (15) 

1 - 9 - 1 - 9 

i{uj s u) = --(^) X(9) sin(- - a), i{d^d) = --(i/)ip) \(9) sin(- + a) 

while the pion condensate, which exists only in the superfluid phase and spontaneously 
breaks the U(l)i symmetry, takes the form, 

z(u l5 d) = ^) Q ^l-\2(9)cos( d -), (fid) = ^Wl-A^sin^) 

Notice that once 9 ^ 0, the P odd condensate iiq^q) appears in addition to the usual P 
even condensate (qq). Similarly, in the superfluid phase at 9 ^ 0, the P even condensate (ud) 
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exists alongside the ordinary P odd pion condensate, i{ifj 5 d). This is a direct consequence 
of explicit parity violation by the 9 term. 

We may also compute the following charge densities from our chiral Lagrangian, 

nj = i$7VV> = /Xl - A(0) 2 ) (16) 

z(u 7 ° 7 5 rf) = -f^i\(6)^l-\i(9)cos(a) (17) 
(u^d) = -ffa\(e)y/l-\*(9)sm(a) (18) 



n A 



All the charge densities vanish in the normal phase. In the superffuid phase, a non-zero 
isospin density appears, rij = -(t^At 3 ^)- This is precisely the density, which one expects to 
induce by applying an isospin chemical potential fij. At 9 = it coincides with the previous 
result I n addition, we also obtain non- vanishing axial charge densities, n A = z(u7°7 5 <i) 



and n A = (u^d). Notice that n A does not vanish already at 9 = 0, nevertheless, it was 
never discussed previously in the literature. The quantity n~ A is the axial charge density, 
corresponding to off-diagonal generators of the SU(2)a group, which is both spontaneously 
and explicitly broken. 

The density n A spontaneously breaks the U(l)j symmetry and, hence, may be considered 
as an order parameter alongside the pion condensate, (ir~) = i(u^ 5 d) Note that there was 
no explicit chemical potential conjugate to n A in the Lagrangian - once U(l)j is already 
spontaneously broken by (tt~), n A is induced automatically. The reader is referred to the 
paper^3[ on N c = Nf = 2 QCD for a few arguments, which intuitively explain why in 
a system with nonvanishing rij and (vr~), the second order parameter n A automatically 
appears. The quantitative behaviour of these two order parameters is somewhat different. 
The pion condensate monotonically increases with \ii after the normal to superffuid phase 
transition, and (n~) — > — |('0V')o f° r ^ On the other hand, the new charge density 
n A first increases after the phase transition, reaches a peak at \ii = ^^m^^ and then 
decreases to for \i ^> m n . Of course, we always consider only fij -C m p . 

We note that the new order parameter n A vanishes, in the limit m q — > 0. We expect that 
in the regime of asymptotically large /xj, where analytical control is present, and both nj 
and (ti~) are believed to be non- vanishing, one can explicitly show that n A will also appear 
once m q ^ is considered. 



IV. THETA DEPENDENCE 

So far we have mostly focused on the \ij dependence at fixed 9. In this section we would 
like to focus more on the 9 dependence, drawing the phase diagram in the (9, fj,j) plane. We 
will also pay particularly careful attention to the physics near 9 = tt. 

We begin by briefly reviewing the well-known 9 dependence at \x = 0. The grand canonical 
potential 0,(9) is, 

n(0 >A i = O) = -flml(9) (19) 
By differentiating Q(9) we can compute correlation functions of GG, 

M = (20) 
d 2 V f 4 , T g 2 GG g 2 GG 
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At fi = we find, 

r 9 2 GG lm u m d . y 

Expressions (122(1 reflect the well-known strong # dependence in the region m u « = m g , 
9 ~ 7r. Let's introduce the asymmetry parameter, e = and assume £ < 1. The 

CP odd order parameter (iGG), though apparently smooth for e ^ 0, experiences a steep 
crossover in the region 1 6* — ?r | ~ e. Correspondingly, the topological susceptibility x has a 
sharp peak around 9 = tt of width A9 ~ e and height x( 7r )/lx(0)| = 1/e. 

Such behaviour of the CP odd order parameter {iGG) strongly suggests that for m u = m d , 
spontaneous breaking of CP symmetry occurs at 9 = tt. This situation, known as Dashen's 
phenomenon, has been extensively studied in QCD with Nj = 3 and Nf = 20,0,0,0,0,0. 
For Nf = 3 with m s ^> m u , m d it is believed that spontaneous CP breaking occurs at 9 = tt 
for \m u - m d \m s < m u m d . 

For Nf = 2, the key observation |4j,|5( is that Dashen's phenomenon is not under complete 
theoretical control in the effective Lagrangian (JTJ). Indeed, for a moment, we fix m u = 
m d . Then, for general 9, the mass term explicitly breaks the symmetry of the effective 
Lagrangian ([!} from SU (2)^ x SU(2)r to SU (2)y. However, for 9 = tt, the mass term in the 
effective Lagrangian vanishes, restoring the symmetry to SU(2)l x SU(2)r and giving rise to 
apparently massless goldstones: m%(9 = tt) = 0. Yet, no such symmetry restoration occurs 
in the fundamental microscopic QCD Lagrangian at 9 = tt. This contradiction is resolved 
by including higher order (quadratic) mass terms in the effective Lagrangian, which would 
explicitly break SU(2)a even at 9 = It is precisely these terms, which control the 

physics of Dashen's phenomenon. 

In this paper we would like to consider two different regimes. In the first regime, one 
may neglect the higher order mass terms by considering fixed 7^ and sufficiently 

small m q . Of course, in such a regime one automatically excludes the regions of parameter 
space where Dashen's transition is realized, and may discuss only the quantitatively steep 
crossover in the normal phase. The second regime that we discuss is obtained by considering 
the exactly degenerate case m u = m d . We show that this second regime exhibits the Dashen's 
transition in the normal phase, which disappears in the superfluid phase. 



A. Crossover Regime 

In this section we discuss the regime in which the leading order chiral Lagrangian (JTJ 
accurately describes the physics for all 9. Here we give only a brief summary of the results 
concerning this regime, for further discussion see |l3| |. 

If the leading order (fT2"|) pion mass at 9 = tt, m%(9 = tt) oc \m u — m d \, is sufficiently 
large one may neglect the higher order mass terms in the effective chiral Lagrangian. For 
any fixed 7^ this is achieved by considering sufficiently small m q . If the higher 

order mass terms are largely saturated by a third quark of mass m Ujd Cm s < ^qcd-, one 
requires, 

\m u - m d \ m u4 m\{9 = 0) 
m u + m d m s 
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FIG. 1: Phase diagram of N f = 2 QCD as a function of m and (9. Here, e = = 0.01. A 

rapid crossover occurs in the normal phase at = tt, which becomes a first order phase transition, 
when m u = nid- 



This condition is, indeed, realized in the true physical world. If, on the other hand, the 
higher order terms are controlled by a light rj' (as motivated by N c — > oo), one needs to 
consider, 

\m u -m d \ m(^) m\{Q = 0) 

m u + m d ^ %M* ~ Ml, 1 ] 

Let us now turn on finite fxj. Once conditions (|23p. are met, all the results of 
previous sections hold for any 9. In particular, the transition to the superfluid phase occurs 
at \i — m 7T (9) (see Fig.HJ). 

Thus, for fx i < m n (6 = 7r) the normal phase is realized for all 9, while for ilj > m n (9 = 0) 
we are entirely in the superfluid phase. Finally, if we fix tii with m- K {9 — 7r) < /// < m 7T (9 = 
0) and vary 9 from to 2-7T we encounter two phase transitions: from normal to superfluid 
phase and then back to normal. Thus, the 9 dependence becomes non-analytic in this region! 
Since the normal to superfluid phase transition is second order, we expect the topological 
susceptibility, \ to be discontinuous across the phase boundary. The transitions between 
normal and superfluid phases occur at 9 = 9 C and 9 = 2n — 9 C , with the critical 9 C given by 
m 7r (0 c ) = hi. 

In the superfluid phase, the free energy density reads, 

m=-\fy r (^^f}- (25) 

Clearly, the 9 dependence in the superfluid phase is different from that in the normal phase 
([19]) . This is most clearly seen by computing topological density and topological suscepti- 
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bility in the superfluid phase, 

g 2 GG m u m d T 2 . 

X = - W^o 2 cos(^). (26) 

The corresponding expressions should be compared with (|2"2"|) describing the normal phase. 
Focusing for a moment on fij > m n (9 = 0), we see that the 9 dependence is very smooth: 
there is no sign of rapid crossover in (iGG) near 9 = it and the large peak in the susceptibility 
X disappears. Moreover, as fij increases, the 9 dependence is suppressed, as expected. This 
smooth 9 dependence at 9 ~ 7r in the superfluid phase should be contrasted with sharp 
behavior in the normal phase discussed above, see eq. ()22j) . As has been explained in detail 
in the parallel study on N c = 2 QCD^I, the disappearance of the "Dashen's crossover" as 
fij increases is accomplished in the following way. First, when m 7r (^ = 7r) < fij < m^ifi = 0), 
the crossover splits into two second order normal to superfluid phase transitions. These phase 
transitions replace the peak in the topological susceptibility x by finite jumps in x at the 
transition points: 

, 2 

sin 2 (9 c ) (27) 



|X(0)| 4/X 

Finally, once /// > m n (6 = 0) no phase transitions can be triggered by varying 9, and the 
"Dashen's crossover" becomes entirely washed out. 

We conclude this section by noting that we can use our results for the topological sus- 
ceptibility x and the chiral condensate (M)ib) to study how the Ward Identities get saturated 
in different phases with arbitrary 9, 0,13 IE 



d ' x{T !^? (0)U " = wj^ M ^ + o(m2) (28) 

0(M 2 ) = ~j d 4 x(T^ 5 M^j(x) tfj l5 M^(0)) conn . 

This Ward Identity is related to the axial anomaly and, thus, should not be affected by 
infra-red effects, such as finite chemical potential. One can explicitly check that our results 
imply that at 9 = 0, the identity (|28j) is, indeed, straightforwardly satisfied both in the 
normal and superfluid phases. However, at 9 ^ 0, in the superfluid phase, one must include 
the 0(M 2 ) term in ()28j) on the same footing as the O(M) term for the Ward Identity to be 
satisfied. The reader is referred to the work where N c = 2 case was discussed in detail. 
In the present case with N c = 3 the saturation of the Ward Identities goes precisely in the 
same way as in [13], and therefore, we do not need to repeat it here. 



B. Phase Transition Regime 

In the present section, we would like to consider the degenerate case m u = rrid = m, which 
has not been discussed in the companion paper This regime is believed to support a 
first order phase transition across 9 = ir at zero chemical potential 0, 0. The discussion of 
the crossover regime in section IVA is highly suggestive of the fact that this phase transition 
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disappears in the superfluid phase. We shall now explicitly demonstrate this claim. We note 
that the point m u = m^, 9 = n might be of importance for lattice fermionsj^, E3) as ft is 
equivalent to a theory where one quark mass is negative and 9 parameter is not explicitly 
present. 

As already mentioned, in QCD with Nf = 2, one needs to include second order mass terms 
in the effective chiral Lagrangian in order to accurately describe physics near m u = rrid, 
9 = it. As argued in [H, the dominant second order mass term is, 



\ 2 (U) = ^/ 7 ^ (ImTr(MU)) 2 (29) 



as it contains 9 dependence different from the leading mass term. Including this term in our 
chiral Lagrangian, we obtain, 

£ = l -flTi{V,UV^) - EReTr(MU) - l 7 ^ (ImTr(MU)) 2 (30) 

Let us review the 9 dependence contained in (|50|l at /// = 0. As is known, the physics 
at hi — in the neighborhood of 9 = tt, crucially depends on the sign of l 7 . A number of 
arguments^ suggest that l 7 is positive, in particular, in the large N c limit, l 7 ~ (HEf. 

n' 

We shall assume l 7 > for the rest of this work. In this case, the static classical minimum 
of (HUH) is given by, U = 1 for < 9 < n and U = -1 for < 9 < 2n. At 9 = it, the classical 
minimum is degenerate: U = ±1, signalling spontaneous breaking of the P, CP symmetries. 
Computing the value of the CP order parameter, (iGG), near 9 = n, 



i- 



g GG m 



327r2 )^ = ±^K^)ol (31) 

This is exactly the result one would derive by naively setting m u = irtd in eq. ()22j) . Thus, 
once m u = rrid, the rapid crossover discussed in the previous section becomes a phase 
transition. 

The 9 dependent mass of the three degenerate goldstones becomes, 

ml(e) = H«M , cos(e/2)| + 2¥n'<M>)l sia , m (32) 

J TV J IT 

We see that the goldstones pick up a small, but non- vanishing, mass at 9 = n, due to the 
Vi term in the chiral Lagrangian. 

Turning on a finite chemical potential, we see that for \fij\ < m ll {9), we are in the 
normal phase, while for \fij\ > m n (9), we are in the superfluid phase (see Fig-EJ) . The static 
minimum U of the Lagrangian ()30|) in both phases is again given by expression (|13|h except 
that now, 

sgn(cos(9/2)) normal phase 
X(9) = { ^(o)cos( g/2 ) superfluid phase (33) 



/x|-m2(7r)sin 2 (6»/2) 



The normal phase at finite \ij again has the same physical properties as at jii = 0. In 
particular, a first order phase transition across 9 = tt persists for \fij\ < m n (9 = ir). 
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FIG. 2: Phase diagram of Nf = 2 QCD for m u = m^. Solid line indicates a first order phase 
transition, while dashed lines indicate second order phase transitions. The region near the triple 
point is subject to further investigation. 



However, once > m n (9 = tt), the Dashen's transition splits into two second order 
normal to superfluid phase transitions. The 9 dependence in the superfluid phase is very 
smooth. In particular, P parity is not spontaneously broken at 9 = tt: one may check, 



(i 



g 2 GG 

327T 2 







(34) 



It is amusing to note, that in the superfluid phase, spontaneous breaking of parity is shifted 
from 9 = 7i to 9 = 0. 

We observe that in the region m n (9 — tt) < Hi < m n {9 = 0), the 9 dependence is 
essentially the same as in a theory with l 7 < and \ij = 0. Indeed, in such a theory, one would 
have instead of a first order phase transition at 9 = tt, two second order phase transitions just 
before and after 9 = tt, accompanied by spontaneous breaking of the SU(2)y symmetry|5J. 
This is similar to the picture that we obtain for l 7 > and m^iji) < fij < m 7r (0), except 
that only the U(l)j subgroup of SU(2)v is broken spontaneously (the other generators of 
SU(2)v are explicitly broken by finite fii). 

Finally, we would like to comment regarding the triple point 9 = tt, = m n (Tr) that 
appears in the phase diagram of Fig. El At this triple point, the set of classical minima 
of the Lagrangian (|30|l presents a sphere S 2 . Such degeneracy is definitely accidental, and 
we expect that it will be lifted by higher order terms in the chiral Lagrangian (most likely 
0(m 3 ), 0(m/x 2 ) terms). These higher order terms also have the potential to change the 
phase diagram in the immediate vicinity of the triple point. However, we believe that once 
we are outside the window, 



10 -d < 



m 



(JLj 



A 



QCD 



m^TT) 



II < 



m 



A 



QCD 



(35) 
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all the results described above are valid. Most importantly, Dashen's phenomenon is present 
in the normal phase and disappears in the superfluid phase. Whether this disappearance 
occurs precisely at the triple point (which would be the most simple scenario) or through 
a more complicated series of phase transitions closely surrounding the triple point is still 
an open question. In order to answer this question one should classify all the terms in the 
effective chiral Lagrangian similar to the classic construction 0]. However, near 9 = ir, the 
dimensional counting rule in such a Lagrangian should be based on the relation m n ~ m q . 
This is in contrast with canonical relation m 2 ~ m q which is the basis for the classification 
scheme presented in We did not attempt to analyze the corresponding problem of 



classification of higher order terms in the effective chiral Lagrangian with fij ^ 0, 9 ^ in 
the present study. As we already stated, outside the region (J35J) the higher order corrections 
can not change our results. 

To conclude the section: we have analyzed the effects of finite fij on the rapid crossover, 
which in the absence of chemical potential occurs at 9 = tt for m u ^ and becomes a 
phase transition once m u = m^. In both cases, the crossover (first order phase transition) 
is replaced by two second order phase transitions as fij increases. We note that if Nf > 2 
for the Dashen's phenomenon to happen one does not require precise equality of the light 
quarks, m u = m^, rather it is sufficient if the quark masses are close enough^, 0]. 2 Our 
remark here is as follows: we believe that in the case Nf > 2 the pattern of the replacement 
of the first order phase transition by two second order phase transitions with increasing \ij 
remains the same as described in the present section. 

Our last remark: Dashen's phenomenon as well as the 9 dependence has been studied 
recently in [19] in a very different approach. The corresponding study had concentrated on 
the weak coupling regime when Euclidean space time volume L is small in comparison with 
the Goldstone mass, L -C m~ . We emphasize that the results presented here are valid in 
the opposite regime L ^> which corresponds to the physically relevant case. 



V. GLUON CONDENSATE 

Having determined the 9 and /// dependence of different condensates and densities con- 
taining the quark degrees of freedom, one can wonder if similar results can be derived 
for the gluon condensate (G 2 ^), which describes the gluon degrees of freedom. As is 
known, the gluon condensate represents the vacuum energy of the ground state in the limit 
m q = 0, fj, = and plays a crucial role in such models as the MIT Bag model, where a 
phenomenological "bag constant" B describes the non- perturbative vacuum energy of the 
system. The question we would ideally want to answer: how will the gluon condensate (G 2 ^) 
(bag constant B) depend on /i, 9 if the system is placed into dense matter? This question 
is relevant for a number of different studies such as the equation of state in the interior of 
neutron stars, see e.g.|20l|. or stability of dense strangeletsjll|. Of course, it is difficult to 
answer this question in QCD at finite baryon density, however, the answer can be easily 
obtained in QCD with fij <C m p , which is the subject of the present work. 



1 This phenomenon, when "naively" higher order corrections in m q start to play a crucial role has been 
previously observed in eq. (|28|l when Ward identities have been analyzed. 

2 It is actually possible that for Nf — 2 the phase transition also occurs already for a very small, but 
non-zero m u — md|(| 
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We work in Minkowski space in this section. We start from the equation for the conformal 
anomaly, 

32tt 2 



% = -^G%G^ + ^Mil; (36) 



where we have taken the standard 1 loop expression for the (3 function and b = yiV c — |iVj = 
Y> f° r N c = 3, Nf = 2. As usual, a perturbative constant is subtracted in expression (|3l)j). 
Now, we can use the effective Lagrangian (JI} to calculate the change in the trace of the 
energy-momentum tensor (O^) due to a finite isospin chemical potential. The energy density 
e and pressure p are obtained from the grand canonical potential Q, 

e = Q + fiini, p = —Q (37) 

Therefore, the conformal anomaly implies, 



b 9 2 na ™n / b 9 2 

v 32tt 2 h ' m ' d W 

(38) 

Here, the subscript on an expectation value means that it is evaluated at fi = m — 0, 9 = 0. 
The good news is that we have already calculated all quantities on the righthand side of 
eq. (J3~8*j) - see expressions (fT5Jl . (fTTJj) . (|23j) . Thus, in the normal phase we obtain, 



(db^G^W - fe^o = -3 m (°)(^)o (39) 



b 9 2 r g ™n _ ,W_ 

'32tt 2 ^ W W 

When 9 = 0, ()39j) reduces to the standard result 0], which was derived in a different manner. 
As expected, {G 2 ) does not depend on /z in the normal phase. The superfluid phase is more 
exciting, 

(^G^Gn,,m,e - (^G^Gno = flA (l + . (40) 

It is instructive to represent the same formula in a somewhat different way, 



\i^^G a tlv G ilua ) ^ mfi - (^^G^G'^ a ) At= o, m ,e - flip] - m l(9)) ~ 2 ^ " 



which makes contact with the fact that in the normal phase, when /z/ < m n {9), the gluon 
condensate does not vary with However, for /zj > m n (9), the dependence of the gluon 
condensate (G 2 ) on /zj in the superfluid phase becomes rather interesting. The condensate 
decreases with /zj for < /z/ < 2 1//4 m 7r and increases afterwards. The qualitative difference 
in the behaviour of the gluon condensate for /zj « and for m n <C /z/ <C m p can be 
explained as follows. Right after the normal to superfluid phase transition occurs, the 
isospin density nj is small and our system can be understood as a weakly interacting gas 
of pions. The pressure of such a gas is negligible compared to the energy density, which 
comes mostly from pion rest mass. Thus, (0^) increases with m and, according to the 
anomaly equation (|36|). (G 2 ) decreases. A similar decrease in (G 2 ^) with baryon density 
is expected to occur in "dilute" nuclear matter (see |22| and review j23|). On the other 
hand, for /zj 3> m n , energy density is approximately equal to pressure, and both are mostly 
due to self-interactions of the pion condensate. Luckily, the effective chiral Lagrangian (JTJ 
gives us control over these self-interactions as long as fij <C m p . Such control is largely 
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absent in corresponding calculations of (G 2 J) in nuclear matter. As Ae ~ Ap, the trace 
(@P decreases and the gluon condensate increases with isospin density. Such behaviour of 
(G 2 ) is quite unusual, as finite quark chemical potentials, on general grounds, are expected 
to suppress the gluons. 



VI. CONCLUSION 

The main purpose of this work was to investigate the phase diagram of Nf = 2 QCD 
at finite 9 parameter and isospin density. We have found that the 9 dependence becomes 
non-analytic: for fixed [ij of order of the pion mass, two phase transitions of the second order 
occur as 9 varies from 2ir to 0. We have also demonstrated the conjecture originally presented 
in |l3| : in the limit of degenerate quark masses, spontaneous P breaking occurs in the normal 
phase, but is absent in the superfluid phase. For m u = m^, a first order transition across 
9 = it is present in the normal phase, but disappears in the superfluid phase by splitting into 
two second order normal to superfluid transitions. The precise details of the neighborhood 
of the triple point where such splitting takes place remain to be determined. 

There are a few more interesting observations which deserve to be mentioned here: 

a) Knowledge of 9 dependence of different condensates allows one to calculate the topological 
susceptibility and other interesting correlation functions as a function of /i. Corresponding 
Ward Identities at nonzero fij are satisfied in a quite nontrivial way, and can be tested on 
the lattice. 

b) Physics of gluon degrees of freedom and fj,j dependence of the gluon condensate can also 
be tested on the lattice. The behavior of the gluon condensate as a function of /i/ is very 
nontrivial, as has been explained in the text. Nevertheless, our prediction is robust in a 
sense that it is based exclusively on the chiral dynamics and no additional assumptions have 
been made to derive the corresponding expression. Our formulae might be useful for the 
lattice simulations when one tries to extrapolate the results to the chiral limit at nonzero 
Hi. 

Finally, we should emphasize that all results presented above are valid only for very 
small chemical potentials fij <C A.qcd when the chiral effective theory is justified. For larger 
chemical potentials we expect a transition to a deconfined phase at fii ~ 5Aqqd HI • 
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